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Abstract: For a class of non-selfadjoint semiclassical operators in dimension one, 
we get a complete asymptotic description of all eigenvalues near a critical value of 
the leading symbol of the operator on the boundary of the pseudospectrum. 

1 Introduction and main result 

In his work pQ, E. B. Davies considers the non-selfadjoint harmonic oscillator, 



By an analytic continuation argument, it was established that the spectrum of P(h) 
consists of the eigenvalues {e iaj/2 h(2k + 1), k = 0, 1, 2, . . .}, while the semiclassical 
pseudospectrum of P(h), defined as the range of the symbol of P(h), p(x,£) = 
£2 _|_ e *a x 2^ on p2 j g e q Ua i t Q sector {z G C;0 < argz < a}. As h — > 0, 

the spectrum accumulates precisely at the origin, which is the only point zq at the 
boundary of the pseudospectrum, for which dp = on p _1 (z ) H R 2 . The angle at 
which the eigenvalues approach the origin is given here by a/2. Now associated with 
the quadratic form p, there is a Hamilton map F, defined by 



with the eigenvalues ±/i, where p, = ie ia ^ 2 . We may therefore reformulate the result 
of PO by saying that the angle at which the spectrum approaches the boundary point 
is given by axg(p/i). In this paper we shall give a straightforward extension of 



USA 



P(h) = h 2 D 2 x + e ta x 2 , 0<a< tt/2. 




this result to what we think is a natural class of non-selfadjoint /i-pseudodifferential 
operators in dimension one. We shall show that, in the semiclassical limit, the di- 
rection of the spectrum near a critical value on the boundary of the pseudospectrum 
is determined by a suitable eigenvalue of the Hamilton map of the Hessian of the 
principal symbol at a corresponding critical point. 

Without any restriction on the dimension, pseudospectra of non-selfadjoint semi- 
classical pseudodifferential operators have been studied in the recent paper which 
extended some of the results of (2|, , [TH]. Following [TH] and [3], let us consider 
the convection-diffusion operator with a quadratic potential, 

P(x, hD x ) = (hD x ) 2 + i(hD x ) + x 2 = e x/2h o UhD x ) 2 + x 2 + ~J o e' x/2h . 

Using the last representation, one can prove that the spectrum of P(x, hD x ) is given 
by {{2k + l)h + |, k = 0, 1, 2 . . .}, while the range of the symbol of P(x, hD x ) fills 
out the region {z G C; Re 2 > (Imz) 2 }. In this case, the spectrum lies strictly inside 
the pseudospectrum, and Theorem 3 in [3| gives general conditions on a point zq at 
the boundary of the pseudospectrum of an operator P = p w (x, hD x ) + 0(h), to be 
away from the spectrum, as h — > 0. The first one is the principal type condition, 

dp ^ along p _1 (z ), (1-1) 

and the second is the following dynamical hypothesis, 

For some A G C, no complete trajectory oiH Re ^x P ) is contained inp _1 (z ). (1.2) 

Asssuming (|1.1|) and (jl.2|) together with a basic analyticity assumption, it is estab- 
lished in |3j that there is a sufficiently small but fixed neighborhood of Zq in C which 
does not intersect the spectrum of P. Back to the one- dimensional case, the purpose 
of this note is to describe the eigenvalue distribution near a boundary point zq, in a 
simplest situation when the condition (jLip fails to hold (so that (jl.2j) fails as well). 
In doing so, as in [3 a , we shall work under the analyticity assumptions, which will 
be essential for our methods. 

Let 

P = P w (x,hD x ;h) 

be the h-Weyl quantization of a symbol P(x, £; h), which is a holomorphic function 
of (x,£) in a tubular neighborhood of R 2 C C 2 , with 

P(x,£;h) = 0(l)m(Re(x,£)) (1.3) 
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there. Here we assume that 1 < m G C°°(R 2 ) is an admissible weight function in 
the sense that for some fixed Co, N > we have 

m(X) <C (X -Y) N m(Y), X,Y eR 2 . 

We next assume that as h — > 0, P(x, £; h) has an expansion in the space of holomor- 
phic functions satisfying the bound (jl.3|) . 

oo 

P(x,Z;h)~'22h j p j (x,Z). 

j=0 

For h > small enough, and when equipped with the domain H(m), the naturally 
defined Sobolev space associated to the weight m, P becomes a closed densely defined 
operator on L 2 (R). 

We shall denote the principal symbol of P by p := po, and make an additional 
assumption of the ellipticity of p near infinity, 

\P(X)\>^^-, \X\>C, |ImX|<± C > 0. (1.4) 
Throughout this paper, we shall work under the assumption that 

m(X) -> oo, \X\ -> oo, X e R 2 . (1.5) 
As before, associated to P, we introduce the semiclassical pseudospectrum 

S(p):=p(R 2 ). (1.6) 

This is a closed set, since p is proper, and we shall assume that E(p) is not all of C. 
It follows then from (|1.5J) that if z is m the complement of then (P — -2o) _1 

exists and is a compact operator on L 2 . The analytic Fredholm theory implies that 
the spectrum of P is discrete and consists of eigenvalues of finite multiplicity. 

Let Zq E dTj{p) and assume that p _1 (z ) H R 2 is a finite collection of points. It 
will in fact be sufficient to treat the case when the pre- image of z is just a single 
point, say (0,0) G R 2 , and in what follows we shall work under this assumption. 
We introduce an exterior cone condition, 

There exist eo > and 6q G R such that (1.7) 
(zq + (0,e )e i{e °- eo ' eo+eo) ) n E(p) = 0. 

Our final assumption is that p — z vanishes at (0, 0) precisely to the second 
order, so that in a neighborhood of this point, we have 

\p(X) - zo\ ~ \X\ 2 . 

The following is the main result of this work. 
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Theorem 1.1 Assume that zq is at the boundary of the semiclassical pseudospec- 
trum of P, defined in (II. and that p~ 1 (z ) R R 2 = (0,0), with p — Zq vanishing 
there precisely to the second order. Furthermore, we make an exterior cone assump- 
tion fll.7|) . Then there exists a G C, \a\ = 1, such that Reap"(0,0) > 0, and the 
spectrum of P in a sufficiently small but fixed neighborhood of z$ in C is given by 

z k = z + G\h{k + ) i ) ] h\+O{h 00 ), keN. 

Here G(q; h) is holomorphic in q G neigh(0, C), and has an asymptotic expansion in 
the space of such functions, as h — > 0, 

oo 

G(q; h)~Y. h ° G M)- 

3=0 

We have G (0) = and 

where fi is the unique eigenvalue of the Hamilton map ofp"(0,0), for which 

Re (a/i/i) > 0. 

We illustrate Theorem 1.1 by comparing the results of a numerical computation 
of small eigenvalues of the operator P = [hD x ) 2 + ex 2 + x 4 , when c = 1 + 3i, with 
the direction of arg(c)/2, given by the quadratic part of the symbol — see Figure 1 
on the next page. 

Acknowledgement. I am grateful to Maciej Zworski for fruitful discussions on the 
subject of the present note. I would also like to thank Caroline Lasser for helpful 
discussions and advice, and Johannes Sjostrand for stimulating comments. The 
partial support of the National Science Foundation under grant DMS-0304970 is 
gratefully acknowledged. 

2 Proof of Theorem 1.1 

Throughout this section, we shall assume, as we may, that zq = 0. Our starting 
point is the following essentially well-known result. 
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Figure 1: Numerical computation of the eigenvalues of P = (hD x ) +cx +x , when 
c = 1 + 3i, for two different values of h. When computing the eigenvalues of P, 
following |14| . we discretized the operator using the Chebyshev spectral method. 
To compare the results of the computation with Theorem 1.1 we have also plotted 
the black solid line of the slope tan(arg(c)/2) 0.7208. 
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Proposition 2.1 Let q(x,£) be a complex-valued quadratic form on R 2 such that 
q(x,£) = precisely when (x,£) = (0,0). Then the range of q on R 2 is either all of 
C or there exists a proper closed convex cone C C which contains the range of q. 



Proof: This result is established in Lemma 3.1 of and we shall recall the proof 
for completeness only. In suitable linear coordinates on R 2 , we may write for some 

q(x, £) = c(x — i£)(x — A£), ImA^O. 
Writing z = x + i£, we find that in the case when ImA < 0, 

q = —z{z + 7z), a 



a c{% — A) 

where 7 = (z + A)(i — A) -1 satisfies I7I < 1. Therefore, 



f^ 2 



Re(ag) = Re (\z\ z + 7 (J) 2 ) = Re ^1 +7^2" J J > °> 

when z 7^ 0, and in this case the range of q is contained in a closed angle of opening 
< 71. In the opposite case, ImA > 0, we find that with some (3 7^ 0, 

(3q = 1 \z\ 2 + (z) 2 , | 7 | < 1. 

It follows that the argument variation of j3q along the circle \z\ — 1 is non- vanishing, 
and hence the range of q on R 2 is all of C. □ 

Remark. According to Theorem 2.1.18 in jB], the range of a general complex- valued 
quadratic form on R 2 , restricted to \z\ = 1, is an ellipse, possibly degenerated to an 
interval or a point. 

Applying Proposition 2.1 to the quadratic form — (p"(0, 0)X, X)/2, which 

begins the Taylor expansion of p at the origin, and using the exterior cone condition 
(jl.7|) . we conclude that since is a boundary point of the pseudospectrum of P, 
the range of P2 is not all of C, for otherwise, the argument variation of p along 
the positively oriented boundary of every small disc centered at 0, would be non- 
vanishing. We then know that for some a G C, \a\ = 1, the real part of ctpi 
is positive definite. After a multiplication of P by a, we may and will assume 
henceforth that a = 1, so that the range of P2 on R 2 is contained in a set of the 
form 

T = {w eC; \lmw\ < CRew}, 

for some C > 0. 
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By polarizaton, P2{X) gives rise to a symmetric bilinear form on C 2 , and then 
we write 

p 2 (X,Y) = a(X,FY), X,YeC 2 , 

where a is the complex symplectic form on C 2 , and F : C 2 — > C 2 is the Hamilton 
map of p 2 — see Section 21.5 of [3] and also jjj, for a systematic discussion of quadratic 
forms on symplectic vector spaces. Since Rep 2 is positive definite on R 2 , it is true 
that F is bijective, and it is well known and easily seen that if /i is an eigenvalue of 
F then so is —fi. When X is an eigenvector of F corresponding to the eigenvalue /x, 
we write 

p 2 (X,X) = a(X,fiX) = 2ifia(ReX,ImX). 

Here the left hand side is equal to p2(ReX, ReX) + p2(lmX, ImX), and then we 
see that either fi/i G T or —fi/i 6 T. In what follows we shall let /i stand for the 
unique eigenvalue of F for which Re (fi/i) > 0. Our goal is to show that as h —>■ 0, 
the spectrum of P near accumulates towards the origin in the direction given by 
the argument of fi/i. 

Remark. In Theorem 3.5 of ^0], the spectrum of the Weyl quantization of a general 
complex valued elliptic quadratic form was computed. When specialized to the case 
at hand, the results of ^0] show that the spectrum of p2(x,hD) is of the form 
{—ifih(2k + 1), k — 0, 1,2,.. .}. Our Theorem 1.1 should therefore be viewed as 
an extension of the one-dimensional version of Theorem 3.5 of jTO] to the case of 
a general analytic /i-pseudodifferential operator. The proof of Theorem 3.5 in [TO] 
employs complex linear canonical transformations, in order to essentially reduce the 
quadratic form globally to a complex multiple of the harmonic oscillator. In our case, 
we shall proceed similarly, but now merely microlocally, working near the critical 
point (0, 0) of the leading symbol of P. 

We shall first use real canonical transformations to simplify p 2 . After a real linear 
symplectic change of coordinates, and a conjugation of P by means of a correspon- 
ding unitary metaplectic operator, we may assume that as (x,£) — > (0,0), 

p(x,0 =p 2 (x,0+O((x,{;f) = ^ + ?)+ l -(ax 2 +be+2cx0 + O((x,0 3 ), A > 0, 

(2-1) 

for some a,b,c G R. Here the quadratic form lmp2 can be diagonalized by means 
of an orthogonal transformation, which may be assumed to be orientation preser- 
ving. Hence this transformation is symplectic, since we are in dimension one, and 
we conclude that after an additional real linear symplectic transformation and a 
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metaplectic conjugation, we may reduce the quadratic part of p to the form in (j2.1j) 
with c = 0. We then write 

P(x, = (* 2 + d?) + 0((x, 3 ), 

where ti = (A + m) _1 (A + zfo), and it is no restriction to assume that Imd > 0. After 
an additional real symplectic dilation in (x, £), the symbol £) takes the form 

P(x,0 = (A + m 2 )M|V2 (a: 2 + e^ 2 ) + O((o;,0 3 ), a = argc? £ [0,tt). (2.2) 

Assuming that a > 0, we shall now look for a small IR-deformation A of R 2 , such 
that the restriction of the quadratic part of p to A will have a constant argument, 
in a neighborhood of (0,0). We first argue formally, and consider Ao = T*(e l6 H), 
where 6 £ (0,7r/4) is to be chosen. Notice that with G(x,£) = 9x£, we can write 
A = exp (iH G )(R 2 ). If £ A then x = e ie y, £ = e~ i6 r] for (y,rj) £ R 2 , and 

using (y,rj) as real symplectic coordinates on A , we get 

1 /2 

p\ Ao = p(exp (iH G )(y, V )) = (A + m 2 )H (e 2 V + e^e^) + 0((y, V ) 3 ). 
We see that if we choose 6 = a/4, then 

p(e W (iH G )(y,ri)) = (A + (y* + ^) + 0((y, V ) 3 ) 

((\ + ia)(\ + ib)) l l 2 , , 9 , 
= U M_ ^ (y 2 + r7 2 ) + 0((y)77)3) 

Here we take the square root with the positive real part. Using the symplectic 
invariance of the Hamilton map F of P2, we then check that here the coefficient 

((A + m)(A + ^)) 1/2 /2 

is equal to fi/i, where \i is the eigenvalue of F, for which Re (/x/i) > 0. 

We shall now introduce a global IR-manifold A which in a very small complex 
neighborhood of (0,0) agrees with T*(e j6l R), and further away from that set it is 
equal to R . The manifold A will also be chosen so that the restriction of p to 
A\{(0,0)} is non- vanishing. We remark that the following construction of A is 
similar to some arguments of jH] and ^2] in the theory of resonances. 

It will be convenient to work with the following FBI-Bargmann transform, 

Tu{x) = Ch- 3 / A J e ilp ^ h u{y) dy, x £ C, C > 0, (2.3) 



S 



where 



<p( x >y) = [y + 1 2~ ^ yx ) 



Notice that lm <p'y > and (p" y 7^ 0, so that <p(x,y) is an admissible phase in 
(|2.3p . Associated to T, there is a complex linear canonical transformation : 
(l/: -<^(x, 2/)) ^ <P*fo y)), given by 

«r(y. ^) = ^= (v - ^ Q/2 r/, -ie- itt/2 2/ + rj) , (y, 77) G C 2 . 

From the general theory ^Tj, we know that kt maps R 2 bijectively onto 

A $0 := I fx, *d x $ (x)J , x G cj , 

where $0 is the strictly subharmonic quadratic form given by 

$0(20 = sup — Im^(x, y). 
yen 

A straightforward computation shows that 

/ s / / x(Rex) 2 1 + sin 2 (a/2) (Imx) 2 . , ; „ NT ^ _ , n iN 

$0 (x) = cos (a/2 + -. \ + sin(a/2)Rea;Imx. (2.4) 

2 cos(a/2) 2 

We remark that $0 is even strictly convex. Another computation shows next that 
the image of T*(e ie R), 6 = a/4, under is of the form A$ 1; where $i( 

When e > 0, we take Xe £ Cq°(C), < x < 1> sucn that % e = 1 when |x| < e 1 
Xe = for I a; I > e 2/,£ , and such that 

d a X e{x) = OJe) \x\~ H , \a\ > 0. 
Then the smooth function 

$ = + (1 - Xe)% 

is strictly convex for e > small enough, and satisfies 

$(x) = <£>i(x), \x\ < e 1/€ , 
= $ (x), tal > e 2/e . 



,1/6 
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In order to decrease the size of the neighborhood of where $ 7^ $oj we introduce 
the strictly convex function 

*i(*)=(^) 2 *(^). ^ :=e2A ' 0<i 7 <l, (2.5) 

so that $,,(x) = $o(^) f° r W > ? 7- Notice that taking 77 sufficiently small, we 
may achieve that ^ — $0 is arbitrarily small in C 1 -norm. We then introduce the 
IR-manifold 

2 

A $r? : f = -d^, 
and remark that along A$ , we have 

{ = {W = x , ^ + (l - X, ?^ + 0(e W), 



T] J i dx \ \ T] J J i dx 

uniformly in r]. We claim now that with p = p o , 

\p(x,£(x)\ > l\x\\ C>0, (2.6) 

if x G neigh(0, C) and e > is small enough. Indeed, using ()2.2j) and computing the 
inverse of Kt, we find that 

p(x, = 2((A + ia)(X + ib)f'Hxi + 0((x, £) 3 ). 

Using the strict convexity of the quadratic forms $0 and $1, we get 



Re(ix£(x)) = X e\x\ 2 + (l~Xe)Re(2xd x %) + 0(e\x\ 2 ) 

= Xe \x\ 2 + (1 - Xe)(x, V$ (x)) + 0(e \x\ 2 ) ~ \x\ 



uniformly in e and 77, and the transversal ellipticity property ()2.6|) follows at once. 
We conclude that for every fixed small rj > 0, the IR-manifold A := (kt) _1 (A$, ; ) is 
equal to R 2 away from an arbitrarily small, previously given neighborhood of (0, 0) 
and it agrees with T*(e* e R) near (0,0). Moreover, it is clear that the restriction of 
p to A vanishes precisely at (0, 0). 

Let us now recall the /i-dependent Hilbert space H(A) ~ L 2 (R), associated with 
the IR-manifold A and equipped with the norm obtained by replacing the weight 
$0 by $ = Q v on the FBI transform side. We also introduce the Sobolev space 
H(A,m) which agrees with H(m) as a set, and whose norm is obtained from the 
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H{m)— norm by the same modification of the weight. Deforming the contour in the 
integral representation of P on the transform side, as explained in [9 , we see that 



P = 0(1) : H(A, m) -> H(A). 



(2.7) 



We summarize the discussion above in the following proposition. 

Proposition 2.2 Let P be as in the introduction and recall that in suitable real 
symplectic coordinates, the leading symbol p of P has the form (12 .2j) . There exists an 
IR-manifold A C C 2 , which coincides with A = T*(e l9 lVj, 8 = a/A, near (0, 0) and 
agrees with R 2 outside an arbitrarily small neighborhood of (0,0), such that p ^ 
along A\{(0,0)}. After applying the canonical transformation k,t, associated with 
the transform (|2.3J) . so that R 2 becomes A$ o; with <£>o defined in (|2.4J) . A becomes A$, 
where $ is a strictly convex smooth function defined in ()2.5j) . For x G neigh(0, C), 
we have 



microlocally unitary near (0,0), such that with P denoting the operator in ()2.7|) . we 
have 



where /i is the eigenvalue of the Hamilton map of the quadratic part of p at (0,0), 
for which Im fi > . 

In what follows we shall be interested in eigenvalues E of P with \E\ = 0(e 2 ), 
when < e 1 is sufficiently small but fixed. The eigenfunctions corresponding to 
such eigenvalues are concentrated in a region where \(x, £)| = 0(e), and therefore 
we introduce a change of variables x = ey. Then 




There exists a Fourier integral operator 



U = e - G{x ' hD * )/h = 0(1) : H{A) -> L 2 (R) 



UP = PU, 

microlocally near (0,0). Here P : £ 2 (R) — > £ 2 (R) has the leading symbol 
p(exp(iH G )(x,0)=p{x,0 = ± (x 2 + e) +O((x,0 3 ), -> (0,0) 



-P(x, hD x ; h) = -P(e(y, hD y ); h), h 
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When viewed as an /i-pseudodifferential operator, e 2 P has a corresponding new 
symbol 

-^P(e(y,v)) ~ -^P( e (y,v)) +hpi(e(y,rj)) + . . . , 

to be considered in a region where \(y, = Oil). 

This scaling reduction together with Proposition 2.2 allows us to reduce the 
further analysis to an analytic /i-pseudodifferential operator P e , microlocally defined 
near (0, 0) G R 2 , whose (not necessarily real) leading symbol has the form 

Pe (x,0 :=po, e (x,Z) = ^±i! + er e (x,0, r e (x,0 = O((x,0 3 )- 
The full symbol is 

oo 

P e (x, 

with Pj, e (%, holomorphic in a fixed complex neighborhood of (x, £) = (0,0). 

From Section 5 of 4\ we recall that for every small enough e, there exists a 
holomorphic canonical transformation 

K e : neigh((0, 0), C 2 ) -> neigh((0, 0), C 2 ), 

such that k £ (0,0) = (0,0), and Im/« e (x, ^) = O(e) when x, £ are real, and such that 

Pe o « e (a:,0 = & 

Here g e {E) is an analytic function of E 1 , with g € (0) = 0, and such that 1m g € (E) = 
0(e) for real E, and Re^(0) > 0. We may also remark that we have in fact 
g € (E) = E + 0(eE). 

Associated to the canonical transformation k £ , we introduce a global IR- manifold 
A e C C 2 , which is e-close to R 2 , equals k € (R 2 ) in a complex neighborhood of (0, 0), 
and agrees with the real phase space further away from this set. Precisely as in jl], 
we then implement K e by means of an elliptic /t-Fourier integral operator U e , 

U e = 0(1) : L 2 (R) - H(A e ). 

Then the action of P e on H(A e ) is, microlocally near (0,0), unitarily equivalent to 
the /i-pseudodifferential operator 

P e = U; l P e U e : L 2 (R) -> L 2 (R), 
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whose complete Weyl symbol has the form 

oo 

P e (x, 

3=0 

with 

Po l6 (ac,0 = 9e 

Continuing to follow |lj (where also further references are given), let us now 
recall how to simplify the lower order terms in an operator whose principal symbol 
is modelled on the one- dimensional harmonic oscillator. In doing so, we consider a 
formal /i-pseudodifferential operator Q(x, hD x ; h), defined microlocally near (0, 0) G 
R 2 , with symbol 

Q(x,& h) ~ q (x,£) + hqi(x,£) + 
where qj are holomorphic in a fixed complex neighborhood of (x, £) = (0,0), and 

with g holomorphic near 0, (?o(0) = 0, g' (0) ^ 0. In Section 5 of |3] we recalled 
that by means of an averaging procedure we can construct 

Mp, & h) ~ a (x, f ) + ha t (x, £) + . . . , 

with all a,j holomorphic in a j-independent neighborhood of (0,0), such that 

e iA ^ hD ^Q(x, hD x ; h ) e - iA ^ hD ^ = g f ^ + fr *? . ^ . (2 . 8 ) 

Here 

oo 

g(E-h)~Y,9AE)h s , 

3=0 

with go as above, and gj(E) holomorphic in a fixed neighborhood of E = 0. As 
in jlj, when deriving ()2.8|) we make use of the fact that if / is a smooth (holo- 
morphic) function on R 2 which depends on (x 2 + £ 2 )/2 only, then f w (x, hD x ) = 
f (\(x 2 + (hD x ) 2 )), for some / with f = f + 0(h) — see [7j for this essentially well- 
known result. 

Applying this discussion to P e , we conclude that there exists 
A e (x, £; h) ~ A , e (x, + hAi, e (x, + . . . , 
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with each A,- e holomorphic in a fixed neighborhood of (0,0) in C 2 , such that after 
a conjugation by e ^(x,hD x -h) ^ ^ p era t or P e becomes 

where 

oo 
3=0 

with all Gj >e holomorphic in q in a j-independent neighborhood of G C. Moreover, 
Go, e = <7e above. 

Combining Proposition 2.2 together with the reductions above, we arrive at the 
following result. 

Proposition 2.3 Let P and p be as in the introduction and recall that we assume 
that the real part of the Hessian of p at (x, £) = (0, 0) is positive definite. Let e > 
be small enough but fixed, and let us view the operator 

^P(x, hD x ) = ^P (e(y,hD y )^j , x = ey, h = 

as an h-pseudodifferential operator, with the leading symbol p e (y,v) = ^P{ e {U-i V))- 
There exists an IR-manifold A C C 2 , with (0, 0) G A ; which coincides with R 2 away 
from a neighborhood o/(0, 0) and is close to R 2 everywhere, such that the restriction 
°fPe to A is elliptic outside an arbitrarily small neighborhood of (0, 0). Furthermore, 
there exists an elliptic Fourier integral operator 

U : if (A) -> L 2 (R), 

such that, microlocally near (0, 0) G A ; 

U (ip) = P.K 

/fere P e = P e (^(l/2)((hD x ) 2 + x 2 );h^ has the complete Weyl symbol, 

oo 

p e (x,e;fc)~x;w'G^ 

3=0 
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with the functions Gj^ e holomorphic in a common neighborhood o/O 6 C. We have 
G , £ (0) = and 

aig (6^(0)) = arg(£)+ 0(e), 

where /i is the unique eigenvalue of the Hamilton map of the Hessian of p at (0,0), 
for which Im \i > 0. 



With Proposition 2.3 available, we are now in a position to solve a suitable Grushin 
problem in the globally defined Hilbert space H(A), and it is clear that the setup 
of the relevant Grushin problem is exactly the same as in 0j. (See also for a 
recent general presentation of Grushin problem techniques in semiclassical analysis 
and spectral theory.) The analysis of the Grushin problem in question is even 
simplified when compared with [I], as we can use ellipticity along A, when away 
from a neighborhood of (0, 0), while near (0, 0), the operator is reduced microlocally 
to a function of the harmonic oscillator. In conclusion, we obtain the following 
result. 



Proposition 2.4 Let e > be small enough but independent of h. The eigenvalues 
of P in an open disc centered at with radius 0(e 2 ) are of the form 



3=0 

with Gj i€ (q) holomorphic in q G neigh(0, C). We have 

GV(g) = ^2(1 + 0(e)). 

% 

Combining Proposition 2.4 together with a scaling argument as in j^J, [T2*j . 
and j3], we obtain Theorem 1.1. 
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